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Using the quantum Monte Carlo Loop algorithm, we calculate the temperature dependence of the 
uniform susceptibility, the specific heat, the correlation length, the generalized staggered suscep- 
tibility and magnetization of a spin-1/2 chain with random antiferromagnetic and ferromagnetic 
couplings, down to very low temperatures. Our data show a consistent scaling behavior in all the 
quantities and support strongly the conjecture drawn from the approximate real-space renormaliza- 
tion group treatment. A statistical analysis scheme is developed which will be useful for the search 
of scaling behavior in numerical and experimental data of random spin chains. 



I. INTRODUCTION 

Over many decades one-dimensional (ID) quantum 
spin systems have attracted much interest. Despite of 
their simplicity they show a wealth of physical proper- 
ties which provide many keys to the understanding of 
phenomena appearing also in other strongly correlated 
systems. Since the discovery of various quasi-lD mate- 
rials the study of ID spin systems, mainly based on the 
Heisenberg model and its variations, is also of experi- 
mental relevance. Examples of such materials include 
the so-called NINO (Ni(C ft Hi N 2 ) 2 NO 2 (Cl) A), NENP 
(Ni(C 2 H 8 N 2 )2N0 2 (Cl)4))Jffland Sr 3 CuPtO 6 .0In partic- 
ular, the latter system belongs to a class of compounds 
which is compositionally very flexible. It is possible to 
produce compounds of the form Sr 3 MN06 in various 
combinations with M=Cu, Mg, Zn, Yb, Na, Ca, Co, and 
N=Pt, Ir, Rh, Bi. Other examples are the spin ladder 
systems, such as Sr n _iCu„ + i0 2 n (n + 1, the number of 
ladder legs) or CaV 2 5 .Q 

Disorder effects play a particularly important role in 
quasi-lD systems, because even small deviations from 
regularity often destabilize the pure phases £1 Real exper- 
imental systems naturally contain impurities and other 
types of disorder. Therefore it is very important to un- 
derstand the influence of disorder on the properties of 
such systems in order to interpret experimental results. 

This is one reason why in recent years random spin sys- 
tems have been investigated intensively. But also from 
the theoretical point of view, random spin systems are in- 
teresting since they provide a simple model to study the 
interplay between quantum effects and disorder. There 
are various realizations of quasi-one-dimensional random 
spin systems in nature. One class is represented by 
compounds like Sr 3 CuPti_ :E Ir ;r 06.B While the pure com- 
pounds Sr 3 CuPt0 6 (x = 0) and Sr 3 CuIr0 6 (x = 1) 
are antiferromagnetic (AFM) and ferromagnetic (FM), 
respectively, the alloy Sr 3 CuPti_.2;Ir x 06 contains both 
AFM and FM couplings whose fraction is simply re- 
lated to the concentration x of Ir. A corresponding 
minimal model of this type is the spin-1/2 Heisenberg 



chain, where the nearest neighbor exchange coupling is 
+J or —J with certain probabilities .0 Another related 
example is realized in the low-temperature regime of ran- 
domly depleted AFM spin-1/2 Heisenberg ladders. Two- 
leg Heisenberg ladders have a resonating valence bond 
ground state which has short-range singlet correlation 
and a spin excitation gapJj If spins are depleted at ran- 
dom the low-energy properties are changed drastically. 
Each site which lost a spin is accompanied by an effective 
spin 1/2, and the residual interaction among these spins 
is randomly FMjQr AFM with a wide distribution of cou- 
pling strengths .Em Recently, it was also shown that ran- 
dom AFM spin-1 chains including next-nearest neighbor 
couplings can generate effective ferromagnetic couplings 
in the low-energy regime sp,that they should behave sim- 
ilarly at low-temperatureJ13 

In this work we focus on random FM-AFM spin chains 
with the full spin rotation symmetry. One of the most 
insightful techniques to study such systems i§-,tjhe real 
space renormalization group (RSRG) method.ErtZI This 
method-was introduced to study random AFM spin-1/2 
chains,tll~tJ and was recently adapted fAi-the study of 
the class of systems introduced above .EjO The basic 
idea of the RSRG method is the iterative decimation of 
degrees of freedom by successively integrating out the 
strongest bonds in the spin chain. In this way the dis- 
tribution of coupling strengths and spin sizes is renor- 
malized. For many cases an universal fixed point distri- 
bution is reached in the low-energy limit. The random 
AFM spin-1/2 chain belongs to a universality class differ- 
ent from that of the random FM-AFM spin chain. The 
former is characterized by a random singlet phase where 
each spin tends to form a singlet with one other spin even 
over very large distances. For the FM-AFM spin chain, 
however, spins correlate to form effective spins whose-av- 
erage size grows with lowering of the energy scale.t£Hla A 
result of the RSRG method is that the average number 
of original spins included in a single effective spin for a 
given (energy or) temperature T scales as. l . f<: .T~ 2a , and 
the average spin size S cx T~" for T — ► 0.I13O 

We start our analysis with a brief review of the RSRG 
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scheme and its resultst§ll3 for the random FM-AFM 
spin chain which shows a universal scaling behavior at 
low-temperatures (section 4.2). The RSRG treatment, 
however, relies on certain approximations which have 
not so far been independently tested. We use, there- 
fore, a different approach to examine the random FM- 
AFM spin-1/2 chain and the validity of— the scaling as- 
sumption underlying the RSRG scheme O For the anal- 
ysis of the numerical data obtained by the continuous 
time versi&n-of the quantum Monte Carlo (QMC) loop 
algorithn£30 we introduce a statistical description of 
the low-temperature properties (section 4.3). Various 
quantities such as the uniform susceptibility, specific 
heat, generalized staggered susceptibility, correlation and 
generalized staggered magnetization are discussed in sec- 
tion 4.4. We find excellent agreement between our anal-, 
ysis and the RSRG treatment by Westerberg et aLEj'EZI 
Beyond the consistent evaluation of the scaling exponent 
a, our simulation gives quantities mentioned above over 
a wide range of temperatures and our discussion provides 
a technique to analyze (numerical or experimental) data 
for the low-temperature scaling regime. 

II. REAL SPACE RENORMALIZATION GROUP 
TREATMENT 

The real space, renprmalization group scheme used by 
Westerberg et a/.E3il3 can be applied to the large class of 
random spin systems with a Hamiltonian of the form 

h = y^^jjSj ■ Sj+i, (i) 

i 

where the strength and sign of Ji as well as the size of the 
spins Si are randomly distributed according to the distri- 
butions P(J) and P(S), respectively. The basic idea of 
the RSRG method is the iterative decimation of degrees 
of freedom by integrating out successively the strongest 
bonds in the spin chain. The RG step will now be de- 
scribed in more detail. 

For this purpose a bond in the chain is defined as two 
neighboring spins and the coupling J connecting them. 
If a bond is isolated from the rest of the chain, it would 
form a local ground state of maximum (J < 0) or min- 
imum (J > 0) spin with an energy gap A to the first 
excited multiplet. For a FM bond A = -J(S L + S R ), 
while for an AFM bond A = J(\Sl - Sr\ + 1), where 
Sl and Sr are the left and right spins of the bond, re- 
spectively. Now we focus on the strongest bond in the 
chain, defined as the bond with the largest gap A — Ao. 
If the distribution of gaps is broad (this is the case if 
the coupling strength and/or the spin-size distribution is 
broad), the gaps of the two neighboring bonds, Ai and 
A2, are typically much smaller than Ao so that the two 
spins Sl and Sr, to a good approximation, lock into their 
local ground state. Consequently, the bond A , Sl,Sr is 



replaced by a single effective spin S' — \Sl ± Sr\ rep- 
resenting the local ground state of minimum (AFM) or 
maximum (FM) spin. The weaker neighboring bonds are 
then taken into account perturbatively in Ai^/Ao, lead- 
ing to an effective interaction between the spins Si, S', 
and £2 (see Fig. |l|). The spin rotation symmetry is pre- 
served in this procedure, and to first order in A1.2/A0 no 
next-nearest neighbor interactions are generated. At this 
point, however, it is important to note that the crucial 
assumption of Ai, A2 -C Ao is probable, but uncon- 
trolled and independent tests are necessary to examine 
whether this assumption and, hence, the RSRG results 
are correct (section 4.3 and 4.4). 

Integrating out the strongest bond successively in the 
manner described above preserves the the form of the 
Hamiltonian Eq. ([!]) but changes the coupling strength 
distribution, the spin distribution and in particular re- 
duces the energy scale Ao- We expect that after many 
iteratively performed decimation steps the distributions 
flow to a fixed-point distributions with rpi'jicrsal scaling 
behavior. According to Westerberg et al.^BxJ these fixed 
point distributions are universal for a large class of ini- 
tial P(J) and P(S), namely if the initial gap distribution 
P(A) is regular or less singular than P(A) ~ A -0 - 7 . 




The RSRG results00 show that_for this whole class 
of distributions the average number I of original spins in- 
cluded in one effective large spin diverges with Aq 2a as 
Ao — > with an universal exponent a = 0.22 ±0.01. The 
total spin of an effective large spin is the sum of the par- 
ticipating original spins, where each spin enters the sum 
with the (opposite) sign as its neighbor if the coupling 
is ferromagnetic (antiferromagnetic). By a random walk 
argument one finds for the average total quantum spin 
number S 

S(xVl(xAQ a . (2) 
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FIG. 2. Possible cluster forming at given temperature T. The cluster lengths l\, h, 
an average value I = X 2 (T / Jo)~ 2a 



are exponentially distributed with 



At finite temperature T the RG flows will stop at A ~ 
ksT. All bonds with a gap larger than Ao are essentially 
frozen to form effective large spins which are, in turn, 
more or less independent from each other. These free ef- 
fective spins lead to a Curie like behavior of the magnetic 
uniform susceptibility \ per unit length for T — ► 0, since 
the scaling behavior of >5 2 and I is identical, 

s__2L£ I. (3) 

L 3k B T I TAq 2q T y ' 

For the entropy a (per original spin) we can argue that 
only the free effective spins contribute such that 

a = | ln(25 + 1) oc T 2 "lnT (4) 

KB 

in leading order, which is identical to that of the specific 
heat C per spin. 

In summary we have the following physical picture of 
the RSRG scheme. At a given finite temperature T, the 
original spins are grouped in completely correlated clus- 
ters forming an effective spin. These clusters are inde- 
pendent and have an average length I oc T~ 2a . This will 
be the basis of our statistical approach for the discussion 
of the thermodynamic properties. 

III. METHODS 

A. Quantum Monte Carlo calculations 

For the numerical simulation we use as model a FM- 
AFM spin- 1/2 chain 

H=J2JiSiSi+i, (5) 

i 

where Si is the i th spin of the chain and Ji, the cou- 
pling strength of the i th bond, takes random values of 
both positive and negative sign. We assume the bond 
distribution 

P(J) = {t° ^ 0<J<J ° (6) 
|^ otherwise, 

where Jo is the maximal coupling setting the energy scale. 



According to Westerberg et the low energy be- 

havior of the random FM-AFM chain is independent of 
the initial distribution as long as J Va P(J) is regular for 
J — > with y c w 0.7. Our choice particular distribution 
[Eq. (^)] is guided by the feasibility for the numerical 
analysis of the scaling regime. This does not restrict 
the validity of the conclusions since the scaling behavior 
should be universal for all distributions in above class. 
Distribution functions closer to those of the real systems 
mentioned in the introduction, on the other hand, while 
interesting ta_study in connection with experiments on 
real material 2 ^, are not suitable for an accurate analy- 
sis of the scaling regime, since the scaling regime occurs 
at much lower temperature relative to the characteristic 
coupling strength Jo. nrn 

Using the QMC loop methodBEl, we simulate 100 
samples of spin chains of 400 sites in a temperature range 
down to temperatures as low as T = Jo/1000. In the 
QMC loop algorithm the calculations are performed di- 
rectly with a Trotter time interval At = 0, so no ex- 
trapolation in the Trotter time interval is necessary (see 
section 2.7). To obtain good statistics, we consider 100 
different realizations of random coupling spin chains, fol- 
lowing the distribution in Eq. (^|). The temperature 
dependences of the physical properties are calculated for 
each chain separately, and then averaged over the 100 
samples. 

B. Statistical cluster analysis 

The fact that the correlation of spins occurs within 
definite clusters permits the use of a special statistical 
approach for the discussion of the thermodynamic prop- 
erties. Original spin-l/2s belonging to the same clus- 
ter are completely correlated, while correlations between 
original spins belonging to different cluster can assumed 
to be zero (c.f. Fig. ||). Therefore, the contributions of 
these frozen clusters to a physical observable A can be 
considered independently and be estimated simply as the 
zero-temperature (ground state) value within each clus- 
ter (depending only on their length I). 

Performing a statistical analysis of the cluster length 
I allows to go beyond the leading order in the tempera- 
ture dependence of \ an d cr, as predicted by the RSRG 
[Eq. (|3|) and Eq. (Q)], as well as to calculate also the T- 
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dependence of other physical quantities. Namely, the low 
temperature dependence of physical observable A per site 
in a random AFM-FM spin-1/2 chain of length L 3> 1 
can be expressed as 

oo 

A(T) = z X> T (0M<, (7) 
1=1 

where nr(l) is the number of clusters of length / and Ai 
is the zero-temperature (ground state) value of A per site 
for such a cluster. As L goes to infinity, we define the 
probability distribution pr(l) of the cluster length I to ap- 
pear at given temperature T, as pt(1) = n T(l)/J2i "r(0_- 
The total number of clusters ^2 1 tit (I) is equal to L/l, 
where I is the average cluster length. Hence 

oo 

A(T) = jJ2pt(1)IAi. (8) 

As mentioned above, the RSRG approach predicts a pow- 
erlaw scaling behavior for the temperature dependence of 
I 



where A is a dimensionless proportionality factor. 

Let us now consider the distribution of cluster lengths I 
for a given temperature. It is natural to assume that the 
bonds freeze uncorrelated to each other, i.e. each bond 
is frozen with a certain probability independent of the 
location within the chain. Thus, the distribution pt(1) of 
the cluster length I has an exponential form. Its average 
value I is given by Eq. (^|) and, hence, we will write 
Pl(l) = pr(l) in the following. The distribution pj(l) has 
to be normalized for a sum over / from 1 to oo and reads 
in its discrete form 

p T {l) = {e K -l)e~ K \ with k = log(l + j—j)- (10) 

The distribution pi(S) of the spin size of a cluster of 
fixed length I also plays an important role, since the 
contribution Ai to many physical observables can be de- 
termined using pi(S). The effective spin of a cluster of 
length I is equal to the ground state spin quantum num- 
ber S. For spin-1/2 degrees of freedom coupled randomly 
by FM or AFM bonds the probability distribution for 
these quantum numbers is given by (an explanation is 
presented below) 

V 1 

Pl{s) = (i/2 + sy.(i/2-sy.¥ {2 ~ Ss ' o) > 

where for even I, S is integer with < S < 1/2 and 
for odd I, S is half- integer with | < S < 1/2. Because 
the system contains no frustrating couplings, the spin 
quantum number S is identical to the total spin of the 
corresponding completely correlated cluster of classical 



spins. The total spin of a classical cluster of length I, 
however, can easily be calculated. It is the sum of all orig- 
inal spins 1/2, where each spin enters the sum with the 
same (opposite) sign as its neighbor if the bond is ferro- 
magnetic (antiferromagnetic). The summation therefore 
can be visualized by a random walk on a ID lattice. A 
ferromagnetic (antiferromagnetic) bond corresponds to a 
step of one lattice constant to the left (right). The ran- 
dom walk contains of / steps and the distance between 
the end point of the random walk and the origin in the 
ID lattice is equal to 2S and the distribution pi(S) [Eq. 
(fil"|)] gives the probability for such a random walk. 

The two distributions pj{l) [Eq. (0)] a,nd_pi(S) [Eq. 
(|TT|)] together with the scaling behavior of I [Eq. (||)] 
are sufficient to calculate the low-temperature thermody- 
namics of various physical quantities. This will be done 
in the next section for the uniform susceptibility, specific 
heat, entropy and correlation length. 

IV. RESULTS 
A. Intermediate temperature regime 

Before turning to the main point of this chapter, the 
study of the universal low-T scaling regime, we will con- 
centrate on the intermediate temperature range of the 
random FM-AFM spin-1/2 chains. The behavior in this 
range and the way entering the universal low-T are both 
strongly dependent on the initial coupling distribution 
function P(J). To study these points, in this subsection, 
we choose different initial coupling distribution functions 
P{J) having a width f^ o dJ|J|P(J) = Jo/2 and being 
normalized to one (J^° dJP(J) = 1). 

First we investigate the uniform Curie constant T%(T), 
as a measure for the still uncorrelated spins at a given 
temperature, considering random FM-AFM spin chains 
of four sites with open boundary conditions. For these 
systems the bond average can be performed exactly and 
the calculations take only little time. Thus a more ex- 
tended study of the dependence of the behavior on the 
initial coupling distribution P(J) in the intermediate 
temperature range is possible. It is clear that finite size 
effects allow only a qualitative study of this problem, 
but the main characteristic features can already be seen 
in the 4-site chain. 

The principal process leading to the dependence of the 
uniform Curie constant on the initial coupling distribu- 
tion function is the following: For a given temperature 
the bonds which had a coupling J > T are frozen either 
into a singlet (AFM bond) or a triplet (FM bond) config- 
uration. Thus, an uncorrelated bond contributes 1/2 to 
T\, a triplet bond 2/3 and a singlet zero (we set the Bohr 
magneton to unity). Thus, the average contribution per 
bond is reduced if the FM and AFM couplings occur with 
the same probability, as considered here. The decrease 
of T\ corresponds to the rate of freezing of bonds as 
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the temperature is lowered, i.e. the constant distribution 
function leads to a linear temperature dependence, as 
observed for the box distribution (see Fig. || and Fig. |). 
For the roof distribution (shown in the inset of Fig. |), 
compared with the box distribution, the deviation from 
the free spin-1/2 limit of 1/4 starts at higher tempera- 
ture, since the largest possible coupling in the chain is 
larger for the roof than for the box distribution. In the 
intermediate temperature range (O.lJo ^ T < 0.5Jo) the 
decrease of the Curie constant for the roof distribution 
with decreasing T is first smaller than for the box distri- 
bution, but then increases and entering the low temper- 
ature regime it is again larger than for the box distribu- 
tion, as expected from the shape of both initial coupling 
distributions (see Fig. |^). Here, the low-T regime means 
the temperature range, where not only single bonds are 
frozen, but all four spins start to be correlated and renor- 
malization effects of the distribution function becomes 
visible. 
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FIG. 3. Temperature dependence of the uniform Curie con- 
stant T\{T) for 4-site random spin chains with initial coupling 
distribution function P(J) shown in the inset (box and roof 
form) . 

In a second step we analyze the effect of an initial 
coupling distribution P(J) which is no longer broad but 
rather peaked at some values J ^ 0. In Fig. [|the temper- 
ature dependence of the Curie constants for 4-site chains 
with initial coupling distributions 



P(J) = { i3X X(\-5)Jo<\J\<{\ + b)Jo 
\ otherwise, 



(12) 



for different 6 = 0, 0.05, 0.25, 0.4, 0.5 are shown. The case 
(5 = corresponds to a random chain with couplings 
±Jo/2, while the case with S = 0.5 represents the box 
distribution Eq. (||). The latter was already discussed 
above. The other broad distribution S = 0.4 leads to a 
very similar behavior in T\ to the box distribution and 



the small differences can be explained by the same argu- 
ment as before. The random spin chains with peaked ini- 
tial P(J) (S — 0.05 and 0), however, behave completely 
different in the intermediate temperature range. As all 
couplings are approximately of the same strength, the 
picture of successively freezing of the strongest bond with 
lowering temperature no longer applies. According to Fu- 
rusaki et alu in such random chains, in the intermediate 
temperature regime all sequences of only AFM and of 
only FM bonds lock together, forming new effective spins 
of random spin size and which are coupled very weakly 
according to a new and broad coupling distribution. This 
forming of weakly interacting effective spins is indicated 
by the plateau in the uniform Curie constant (see Fig. ||). 
As the temperature is lowered further, the effective spins 
start to correlate and the Curie constant decreases again. 
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FIG. 4. Temperature dependence of the uniform Curie con- 
stant T\{T) for 4-site random spin chains with initial coupling 
distribution function Eq. (|l^ ) shown in the inset. The dif- 
ference between the Curie constant for the case S = 0.05 and 
5 = are smaller then the symbols. 

The study of the 4-site chains, giving an idea of how 
the universal scaling regime is reached for different P(J), 
clearly shows the advantages of our choice of a broad ini- 
tial distribution P(J) [Eq. @] for the numerical QMC 
study of the low-temperature scaling regime. Reducing 
the energy scale, the broad distributions are renormal- 
ized quite rapidly to the fixed point coupling distribu- 
tion, while in the case of the peaked distributions, first 
new effective spins have to be formed before entering the 
universal low-temperature scaling regime. 

Whereas the thermodynamic properties in the interme- 
diate temperature range and the crossover temperature 
to the scaling regime depend strongly on the initial cou- 
pling distribution, the low-temperature regime of random 
FM-AFM spin-1/2 chains itself should show universal be- 
havior. This will be analyzed in the following, based on 
the scaling assumption of the RSRG method, using the 
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statistical cluster analysis, introduced in section 4.3, as 
well as the QMC results for the box distribution Eq. @ . 

B. Low-temperature Curie behavior 

An important result of the RSRG approach is that the 
uniform susceptibility approaches a Curie-like behavior 
in the limit T — > 0. It is possible to determine the lim- 
iting Curie constant from the initial spin-size distribu- 
tion and the ratio of FM to AFM bondsE30. The low- 
temperature deviations from this Curie constant give im- 
portant information on the scaling properties mentioned 
above. For the scaling analysis we calculate the prod- 
uct T\(T) with the QMC simulation over a temperature 
range, Jo/1000 < T < Jo as shown in Fig. ||. According 
to the RSRG scheme, with lowering temperature gradu- 
ally the spins correlate within clusters of growing length 
I (see Fig. |J). Each such cluster represents an effective 
spin. In our case this leads to an effective diminishing 
of the spin degrees of freedom such that Tx(T), as a 
measure for the still uncorrelated spins at a given tem- 
perature, decreases monotonically from the large- T value 
of 1/4 (the free spin- 1/2 limit). The linear dependence 
in the intermediate temperature range reflects the initial 
bond distribution P(J) in Eq. (jfy, as discussed above for 
the 4-site chains. 
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FIG. 5. QMC results for the temperature dependence of 
the uniform Curie constant of random FM-AFM spin- 1/2 
chains of length L — 400 with a bond distribution according 
to (^). Also shown are the results for L = 2, 3 and 4 cal- 
culated exactly by integrating over the distribution for each 
bond. For the small clusters open boundary conditions are 
used to avoid strong frustrations. 

We now turn to the analysis of the low-temperature be- 
havior of x- Since each cluster behaves independently we 
can calculate its contribution to the susceptibility and av- 
erage over all cluster sizes and spin configurations. Thus 



for clusters of the average length / the Curie "constant" 
is given by 

T X (T) = 7 J2pi(1)Ci, (13) 
L i=i 

where Ci is the average value of Curie constant per site for 
a finite, completely correlated cluster of length I which 
is given by the average over all possible spin quantum 
numbers appearing in such a cluster 

i/a 

Ci = -J2M S ) S ( S + 1 )- ( 14 ) 

If I ^> 1, it is justified that we use the distribution func- 
tions pi(S) and Pi(l) in their continuum approximations 

pi{S ) = 2J^e-^ 

(15) 

pi(i) = ^e-mi-D 

for I > 1, as approximation of Eq. ([H^) and Eq. ( pd| ) 
- a comparison, done in Appendix A, of the following 
results with numerical results obtained using the discrete 
distribution shows only minor deviations even for small 
I. These two distributions [Eq. (|l5|)1 are now used to 
calculate the low-temperature behavior of Tx 



T X {T) = l J™ dlp T (l)lCi 
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FIG. 6. Uniform Curie constant per spin as a function of 
T a for different a. The best linear behavior [Eq. (|l^)] is 
found for a = 0.21 ± 0.02. The solid line shows a fit of form 
( p^| ) for a — 0.21 and fixed intercept 1/12. This leads to a 
value of A = 1.0 ± 0.1 . 
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Fitting the two leading terms of Eq. ( |l6| ) to the QMC 
data gives an estimate of the scaling exponent a and the 
proportionality factor A. Fig. |^ shows T\ as a function 
of T a for different a. Using the knowledge of the exact 
zero-temperature value 1/12 of T\ we obtain from this 
plot a = 0.21 ± 0.02--in good agreement with the RSRG 
result (0.22 ± 0.01).IIZl The proportionality factor is de- 
termined as A = 1.0 ± 0.1. These results confirm that 
our QMC results for chains of length L = 400 are not af- 
fected by finite-size effects, because according to Eq. @ 
the average length of the correlated clusters, even at the 
lowest simulated temperature Jo/ 1000, is only I 18, 
and therefore much smaller than L. 



C. Specific heat and entropy 

The specific heat C is determined as the numerical 
derivative of the internal energy calculated in the QMC 
simulations. The errors in the specific heat increase sig- 
nificantly with lowering temperature. Nevertheless, C 
could be determined reliably down to T « J o /150. The 
QMC results for the specific heat and the ratio C/T as 
a function of T are shown in Fig. 
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FIG. 7. Specific heat divided by temperature T as a func- 
tion of T. The circles are the QMC results, while the solid 
line was calculated as the derivative of Eq. (hq) with respect 
to T. The inset shows the QMC results for the specific heat 
as a function of temperature over a large temperature range. 
The error bars are smaller than the symbols. 

Since C /T is the derivative of the entropy a with re- 
spect to T, we may compare it to a calculation of a by 
our statistical cluster analysis. Again we assume that for 
given temperature each correlated cluster behaves inde- 
pendently Hence a cluster of length I contributes 



<7, = £]/>,(S)log(2,S+l). 



(17) 



to the entropy. Here it is necessary to use the discrete ex- 
pressions for the distribution functions because large de- 
viations arise in the continuum approach using Eq. (|l5|), 
instead of Eq. ( |l0| ) and Eq. (11) (see Appendix A). 
The entropy per site of the infinite chain is 
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(18) 



i=i 



where I is given by Eq. (||) . Using the values of a and 
A determined above, the entropy and hence C/T may be 
calculated from Eq. ( |l8| ) . The result is shown as the solid 
line in Fig. fj]. We find good agreement with the QMC re- 
sults for the low-temperature range. We emphasize that 
this is no fit since all free parameters have been deter- 
mined via the uniform susceptibility. At higher temper- 
atures this approach fails, because the scaling behavior 
Eq.([)J) of r is no longer valid. 

As a further consistency check we estimate the area 
below the curve C/T from T = to oo using a (oo) — 
cr(T*) + f™dT(C/T), where we set T* = 0.03J - The 
first term is obtained using Eq. ([l8]) and the integral 
from T* to oo is determined numerically, using the QMC 
results and a high temperature expansion up to third or- 
der. We find cr(oo) = 0.68±0.01, in very good agreement 
with the expected result In 2. The fraction a(T*) of the 
entropy is large, contributing approximately 25% of the 
total. Therefore the fact that we find the correct value 
for q-(oo) is a further convincing test for the correctness 
of Eq. (118) , and of our statistical treatment of the scaling 
regime. 



D. Generalized staggered susceptibility 

The generalized staggered susceptibility is defined here 
as the linear response of the spin chain to a staggered field 
H st whose sign on site j is given by 



I J sgn(-J m ) 



(19) 



in analogy to the regular AFM staggering. Based on 
the discussion of the non-linear sigma model version of 
the random spin chain Nagaosa and coworkers suggested 
that each cluster forms a (classical)-, staggered spin s st 
proportional to the cluster length lf5 



Sst = Ql 



Q = constant 



(20) 



Then the staggered field H st couples to s s t analogous to 
the case of a uniform field and the susceptibility should 
be given by the following low-temperature scaling form, 



7 



1 

Xst oc cx - cx T 



(21) 



Using the proposed Eq. ( |20| ) and performing a similar 
analysis as for the uniform susceptibility, we obtain the 
following form 



Xst 



A j.-(l+2a) , 

3 



C-2C 2 



0{T 



-l+2a 



)• (22) 



However, no £ can be found such that Eq. ( p2| ) fits our 
QMC data. Indeed, the low-temperature behavior of Xst 
is governed by a leading powerlaw Xst T^ 1 with 7 = 
1.17±0.01 which is clearly different from the anticipated 
1 + 2a. One may argue that the temperatures covered 
by the QMC simulations are not low enough and that 
only at lower temperature the exponent will approach 
the value 1 + 2a. The numerical data, however, show no 
indication for such an increase of the exponent down to 
temperatures as low as T — Jo/1000. 
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FIG. 8. Double logarithmic plot of the generalized stag- 
gered susceptibility as a function of temperature. The error 
bars are smaller than the symbols. The dashed line shows a 
divergence with T~ 



-(l+2a) 



This leads to the conclusion, that the assumption in 
Eq. (|2(]) and Eq. (|2l]) is not appropriate, and a revision 
of the ideas on the generalized staggered susceptibility is 
necessary. 

In the following we will follow an alternative path 
which includes quantum effects ignored above and finally 
gives an overall consistent view of this problem. We con- 
sider the application of the staggered field H st as a per- 
turbation described by addition the following term to the 
Hamiltonian, 



V^t — —HstSst — 



(23) 



We now focus on the contribution of a given correlated 
cluster of length I to the generalized staggered suscepti- 
bility. Since the cluster is correlated it is described by its 



ground state I'l'o^ M)) where (S, M) are the spin quan- 
tum numbers (effective spin degree of freedom). This 
state is, in general, not an eigenstate of V^t an d we find 
that (* (S,M)|S st |*oCS , ,M)) grows with cluster size at 
most as Z 1 / 2 (see Appendix B). As we will see below this 
leads to a Curie or sub-Curie law of the generalized stag- 
gered susceptibility. Therefore, we include the second 
order in perturbation, 



AE M = - 



H st (V (S, M)\S«\# (S,M)) 

— \(<f (S,M)\S st \tp n )\ 2 



H' 



eo 



-H st BMl — H st eM2, 



(24) 
(25) 



where \ip n ) denotes the excited state with energy e„ and 
cq is the ground state energy. 

The temperature dependence of the generalized stag- 
gered susceptibility is easily calculated from the free en- 
ergy of the cluster, 



F s = -k B T\nZs, 



Z s = 



s 

E • 

M=-S 



-AE M /k B T 



(26) 



Note that other contributions to the partition function 
Z are of the order e" A "' lfiT and can be ignored (Ao (~ 
fcflT) the excitation gap of the cluster). The generalized 
staggered susceptibility per spin in the cluster is then 
defined as 



Xst 



1 d 2 F 



fls, 



(25- 



1^1 

; M=-S 



C M1 



2e M 2) 



(27) 



Now let us give upper bounds to the two terms. For ejvri 
we find 



Ej/i < 9il p/2 



(28) 



independent of M with a constant g\ and the exponent 
p/2 (see Appendix B). This upper bound is connected 
with the following estimate of ejv/2 which contains two 
elements the denominator and the numerator. The de- 
nominator has a simple lower bound e n — eo > Ao. An 
upper bound of the numerator is obtained if the sum over 
n is extended to run also over the states | v I'o(5', M)), lead- 
ing to 



ei/2 < 



(*o(S,M)\S st S st \y Q (S,M)) 
Ao 



J- ]T TiTjiMS, M)\SiS*\* (S, M)>. 



(29) 



This leads us to the discussion of the generalized 
staggered correlation function which we assume to be 
bounded as 



8 



TiTjiMS, M)\SfS*\MS, AO) < g 2 l f 



(30) 



hi 



where 172 is a constant and p a scaling exponent. It is 
easy to see that p = 2 for "long-range ordered" clusters 
and p = 1 for short range correlation. An intermediate 
p indicates a correlation function which decays with a 
powerlaw, 



Io(r) = 71T <+r (* (5, M)|S? 5* +r |* (5, M)) 



(31) 



for large r and i what leads to p = 2—7/. Combining these 
estimates the generalized staggered susceptibility gets the 
following upper bound from the two terms, Xet = Xi +X2, 



Xi 



E 

M 



-Ml 



(25+ l)ZfcflT 



. < , )il3 t^ocT-( 1 + 2Q( "- 1 )), 



X2 = E 



6M2 



,1/ 



(25+1)/ 



< const • 



A 



T 



x T -(l+2a(p-l)) 



where we took Aq ~ fcsT and I ~ Z 



-2a 



(32) 
Since 

7] = is a lower bound for 77 we find that the exponent 
of Xst ^ T~ 7 lies between 1 and 1 + 2a. Furthermore 
the knowledge of 7 from our numerical calculation allows 
us now to give an upper bound on the exponent rj (7 = 
1.17+0.01 from Fig. §), 



7 < 1 + 2a(p- 1) 

p = 2 - n 



77 < 1 ^— ^ * 0.(i2 ± 0.02 



2a: 



(33) 



This result is consistent with a recent calculation of the 
correlation function of large clusters using the density 
matrix renarmalization group method by Hikihara and 
coworkersEj Determining the correlation function of the 
ground state averaged over many samples they found a 
powerlaw with 77 ss 0.4 — 0.5. We will also see that our 
following data fit well into this interpretation. Hence 
we can conclude that the correlation is longer ranged 
than in a regular AFM spin- 1/2 chain where the exponent 
is 1. On the other hand, the exponent of Xsti smaller 
than 1 + 2a, is inconsistent with the assumption that 
the ground state of the random FM-AFM chain has long 
range order. 



E. Correlations and generalized staggered 
magnetization 

Finally, we calculate the correlation function 
, L /i+r-l \ 

r M = (75 E S i S i+r II "ga(- Jm) ) ( 34 ) 
i— 1 \ m=i / 



for temperatures Jo/400 < T < Jq, using the QMC al- 
gorithm. For large r and at fixed temperature T, the 
correlation function T(r, T) is found to be rather well de- 
scribed by a pure exponential form 



T(r,T) =R(T)e~ r /^ T \ r » £(T). 



(35) 



£(T) is the correlation length and the prefactor R(T) is r 
independent. Fitting Eq. ( |35| ) to the QMC data, leads to 
an estimate of £(T) and R(T). The correlation length as 
function of temperature is shown in Fig. ^| which diverges 
for T — > with the approximate powerlaw, 



£(T) cx T 



-0.46 



(36) 




FIG. 9. Double logarithmic plot of the correlation length 
£ as function of temperature T. £(T) diverges approximately 
with T -0 ' 46 (dashed line) for T -> 0. The error bars are 
smaller or of order of the symbols. 

Following the concept of our statistical cluster analysis, 
the correlation length £(T), a measure of distance within 
spins are correlated in the random spin chain, should 
be proportional to the average cluster length I. From 
this point of view we therefore expect a low temperature 
behavior according to 



£ oc I cx T 



-2a 



(37) 



Putting in the value a = 0.21 ± 0.02, determined be- 
fore, a divergence of £(T) oc j 1 -o.42±o.o4 j g p rec [i CIl ed, in 
reasonable agreement with the observed divergence, Eq. 
( p6| ) . After considering the length scale of the correlation 
function we now turn to the amplitudes. 

We propose the following scaling behavior for the cor- 
relation function of random FM-AFM spin chains (for 
T — > 0): 



r(r)=ff(r/0, (|»1) 



(38) 







where T(x) is a universal (temperature independent) 
function. 

As a consequence, comparing Eq. ( |35| ) to Eq. (|38|), at 
low temperatures, the prefactor R(T) should behave as 
£(T) V , leading to 

In R(T) - bxR(T ) = z/ln£ - z/lnf 

= -2z/amT + 2z/alnT , (39) 

where R(Tq) and £o is the prefactor of the correlation 
and the correlation length, respectively, at fixed temper- 
ature T [c.f. Eq. (H)]. In Fig. [k] the logarithm of the 
prefactor R is shown as function of In T, at very low tem- 
peratures. The slope — 2aiy is 0.28 ± 0.02, leading to a 
first estimate of v = —0.61 ± 0.06. 




FIG. 10. Double logarithmic plot of the prefactor of the 
correlation function R as function of temperature T [cf. Eq. 
(|f) and Eq. (|§)]. The slope is ~2ua = 0.28 ± 0.02 (dashed 
line), leading to an estimate of v = —0.61 ± 0.06. 

Another possibility to estimate the exponent v is to 
consider the low-temperature behavior of the square of 
the generalized staggered magnetization defined by 



Ml 




2 f°° 

- J drT{r) (L»£»l), 



(40) 

(41) 
(42) 



where we replaced the sum by an integral and assumed 
£ ^> 1 (i.e. very low T). This allows us to neglect the 
term T(0)/2L. 

Using the scaling form Eq. (|3^ ) for T(r) in Eq. (|42| ) 
we obtain, 



Ml = =r r / drr(r/0 



L r +1 y o dxT( X ). 



(43) 
(44) 



As the function T(x) is temperature independent, the 
square of the generalized staggered magnetization should 
diverge as 



Ml cx T +1 cc T 



-2a(V+l) 



(45) 



The QMC data of the square of the generalized stag- 
gered susceptibility is shown in Fig. O. The numerical 
results show a divergence of Ml cx T^- 180 ^ - 002 , From 
this we find for the exponent v « —0.61 which is in good 
agreement with the previous result, deduced from the 
prefactor of the correlation. Note, however, that we esti- 
mate a rather large error of about 15% for z;,as a conse- 
quence of the uncontrolled approximation performed by 
replacing the sum by an integral in Eq. ( fll| ) and neglect- 
ing the term r(0)/2£. 
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FIG. 11. Double logarithmic plot of the square of the gen- 
eralized staggered magnetization Ml as function of temper- 
ature T. Ml diverges with r~° 180±0 ' 002 (dashed line) for 
T ^ 0. 

Our QMC results show that the scaling behavior of the 
correlation function T(r) of the FM AFM random spin 
chains is of the form of Eq. ([58|). The exponent v is 
roughly estimated v = —0.61 ± 0.06 and is a third expo- 
nent (apart from a and 7) describing the low temperature 
scaling regime of random FM-AFM spin chains. 

The exponent v is related to the correlation function 
used for the estimate of the generalized staggered sus- 
ceptibility. We would like to give here a brief argument 
for the scaling form in Eq. (|38|). We consider a given 
distribution of clusters for fixed temperature, and aver- 
age cluster length I with the distribution function in Eq. 
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(|l5|). From the assumption that the clusters are com- 
pletely uncorrelated we derive the following approxima- 
tion for the generalized staggered correlation function for 
distance r ^> I 

poo 

r(r) oc / dlpiiiy-* 1 oc r- rt eT rl1 . (46) 

J r 

Here only clusters of length larger than r can contribute 
and give on the average a value scaling with r as r _7) 
[ground state correlation within the cluster, see Eq. ([H])] . 
Naively we might assume that r\ and v are identical. 
However, the numerical calculations show that the co- 
herence length is shorter than the average cluster length 
(£ k, 21/3) although their temperature dependence is ap- 
parently the same. The form given in Eq. (^) may be 
interpreted as an upper bound which does not account 
properly for the thermal fluctuation effects in the (most 
important) long clusters whose excitation gap is of the 
same order as the temperature. 

V. CONCLUSIONS AND SUMMARY 

In the study of the random FM-AFM spin-1/2 chain, 
the continuous time version of the powerful quantum 
Monte Carlo loop algorithm allowed us to reach temper- 
atures where clear characteristics of the universal low- 
energy properties of this system are observable. The 
characteristic temperature where the scaling behavior be- 
gins depends on the initial distribution of bond strengths. 
The choice of a non-singular distribution has the advan- 
tage that this temperature regime is well accessible with 
our simulation techniques. The key feature for analyz- 
ing the low-temperature data lies in the fact that the 
spins are correlated within clusters whose length grows 
with decreasing temperature. The scaling of the average 
length I — X 2 (T / Jo)~ 2a and the statistics of the effective 
spin distribution allow a very good analysis of the numer- 
ical data. We determine both the prefactor and the ex- 
ponent. The latter is assumed to be universal and agrees 
surprisingly well with the exponent found by the RSRG, 
a = 0.21±0.02E3O These two parameters (a and A) are 
sufficient to describe the uniform susceptibility, entropy 
and specific heat in the whole low-temperature scaling 
regime. 

A property not accessible by the RSRG method men- 
tioned above is the generalized staggered susceptibility. 
In our simulation we show that a new exponent 7 ap- 
pears, describing the low-temperature behavior, Xst oc 
T -7 with a value 7 ~ 1.17. A further exponent v 
arises from the scaling form of the correlation function 
Eq. (|38|). We have shown that v (k, 0.61) is con- 
nected with the ground state correlation function which 
probably exhibits a powerlaw decay over long distances, 
ro(r) oc r~ ri < r~ u . Both exponents are assumed to be 
universal like a, since they are derived from the proper- 
ties in the low-temperature scaling regime. At present, 



however, we cannot estimate the accuracy of the derived 
values. There is also no obvious relation between 7, v 
and a. 

Since in the low-temperature regime the system con- 
sists of effective spin degrees of freedom which become 
more and more classical in their natureo, there is a clear 
tendency towards order in the ground state. The ques- 
tion is whether the remaining quantum effects (fluctu- 
ations) are sufficient to destroy long range order, here. 
Therefore, the important result, besides confirming the 
physical picture obtained by the RSRG scheme, is the 
fact that we do not find any indication that the ground 
state has long range order. The generalized staggered 
susceptibility as well as M s 2 t show too increase too slowly 
at low temperatures. Nevertheless, the correlations must 
be clearly longer ranged than in the uniform AFM quan- 
tum spin-1/2 chain whose ground state staggered cor- 
relation function behaves as To(r) oc r . This finding 
is consistent with a recent density matrix renormaliza- 
tion group study and an improved version of the real 
space renormalization group scheme by Hikihara et aL.tH 
Their treatment, however, suggests that the correlation 
function might reach its real long-range behavior only at 
very long distances which we do not access in our finite 
temperature calculation (I < 20) . 

Among the (unfrustrated) random spin chains with 
complete spin rotation symmetry there are only two dis- 
tinct classes, the ones with the random singlet phase 
fixed point and the others which belong to the class 
studied here.o For the former it was shown by Fisher 
that the RSRG scheme describes the fixed point behav- 
ior exactly.Ej This was not possible so far for the random 
FM-AFM spin chain. Our numerical and statistical anal- 
ysis, however, demonstrates very convincingly the valid- 
ity of the universal scaling assumption. Finally, we would 
like to emphasize that our statistical fitting procedure is 
very suitable and useful to analyze not only numerical, 
but also experimental results of this class of random spin 
systems. 
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APPENDIX A: CONTINUOUS VS. DISCRETE 
DISTRIBUTIONS 

There is one point to which attention has to be paid 
when using the statistical cluster analysis. For each con- 
sidered observable it is necessary to investigate if for its 
calculation the exact form of the distributions [Eq. ( fTo| ) 
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and flllD] have to be taken or if the approximate contin- 
uous forms Eq. (|l^) and the replacement of the sums by 
integrals is also reliable. It is clear that the second case 
is the more favorable one, since it gives simple results for 
the low-temperature behavior of the corresponding phys- 
ical observable. Sums over the discrete distributions, on 
the other hand, usually can only be evaluated numeri- 
cally. Therefore, the continuous forms Eq. (|l^) are pref- 
ered as long as they do not lead to a loss of precision. 
This is shown here for the example of the uniform Curie 
constant and the entropy. 
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FIG. 12. Uniform Curie constant Ci per site of completely 
correlated clusters of length /, calculated exactly (solid points) 
and with large-/ approximation (solid line). The inset shows 
the Curie constant of an infinite chain as function of the av- 
erage cluster length I, calculated exactly (dashed line) and 
with approximate distribution pf and pi(S) (solid line). For 
details, see text. 



We first concentrate on the uniform Curie constant. 
The Curie constant for a completely correlated cluster 
as function of its length is shown in Fig. [l2|. It is first 
calculated exactly (solid points), using Eq. (|14[) and the 
discrete form Eq. (11) of pi(S). These exact results are 
compared to the large-/ approximation of Ci (solid line), 
which is obtained using also Eq. (0), but using the con- 
tinuous approximation Eq. (p"5|). It can be seen that this 
approximation is very good, also for small /. As a con- 
sequence, in Eq. (^|), the approximate large-/ value for 
Ci can be used without loss of precision. The effect of 
replacing the exact discrete distribution pi(l) [Eq. (10)] 
by the continuous one [Eq. (|l5|)] and the sum by an in- 
tegral in Eq. (|l4|), is also investigated. In the inset of 
Fig. [l^ we show the uniform Curie constant of a random 
spin chain as a function of the average cluster length /. 
The dashed line represents the result, obtained by us- 
ing both exact (discrete) distributions [Eq. ( |ll| ) and Eq. 
(|l0|)] and the solid line shows the result, obtained by us- 
ing both the continuous distributions Eq. (|l5|). We find 



excellent agreement for / > 3, what corresponds to tem- 
peratures T < Jo/10. This temperature range covers the 
region, where the cluster distribution can be assumed to 
be exponential and our statistical cluster analysis results 
apply. 

For the entropy, on the other hand, the large-/ approx- 
imation of ai [i.e. using the large-/ approximation Eq . 
( |l5| ) for pi(S) instead of the correct distribution Eq. jTlf ) 
in Eq. (|l8|)] is very bad for small / (see Fig. |l3|) and it 
is not appropriate to use this approximation for the cal- 
culation of a (I). This can be seen in the inset of Fig. |l3| , 
where we plotted the exact value a(l) (dashed line) to- 
gether with the approximated value of <r(Z), calculated 
using the large-/ approximation for ai (dashed line). The 
a(l) based on the continuum approximation agrees with 
the exact one only for / > 15, which corresponds to tem- 
peratures T < Jo/500. In order to compare the statistical 
cluster analysis results with the QMC data, it is there- 
fore necessary to calculate the entropy exactly and not 
to use the large-/ approximation. 
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FIG. 13. Entropy contribution oi of completely correlated 
clusters of length I, calculated exactly (solid points) and with 
large-/ approximation (solid line). The inset shows the en- 
tropy per site of an infinite chain as function of the average 
cluster length /, calculated exactly (dashed line) and with ap- 
proximate distribution p\ and pi(S) (solid line). For details, 
see text. 



APPENDIX B: BOUNDS FOR THE STAGGERED 
MOMENTS 

The aim of this section is to establish the upper bounds 
for ejv/i and &M2 in Eq. (^8|) and (|30|), respectively, for 
a given cluster with the ground state \^q(S, M)) of de- 
generacy 25+1. The generalized staggered correlation 
function T(r) given in Eq. (|3l| ) is assumed to scale with 
the powerlaw r~ n . Therefore we find 
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(MS,M)\S st S at \MS,M)) 

= Eh,- TiTjiMS, M)\S?S?\MS, M)) (47) 

= /„' dr f Q dr'T 0M (\r - r'\) oc Z 2 -' = 

which defines the exponent p. The proportionality factor 
Q2 in Eq. (|3(]) was chosen large enough such that gil p 
acts as an upper bound for all spin quantum numbers M 
for given S and I. 

Now we turn to ejv/i for which we obtain an upper 
bound using Eq. ([47|). 

e 2 n = \(^ (S,M)\S st \^ (S,M))\ 2 

< E„(*o(5,M)|5 8 t|*„)(*„|5 8t |*o(S, M)) (48) 
= (* (S,M)\S st S st \*o(S,M)) < g 2 lP 



where the sum runs over the complete set of basis states 
\ty n ) including, of course, \^o(S,M)}. Taking the square 
root we obtain the upper bound for ejui with g\ > J~g~2- 
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